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Abstract 

In the DGP model, the "self-accelerating" solution is plagued by a ghost instability, which 
makes the solution untenable. This fact as well as all interesting departures from GR are 
fully captured by a four-dimensional effective Lagrangian, valid at distances smaller than 
the present Hubble scale. The 4D effective theory involves a relativistic scalar tt, universally 
coupled to matter and with peculiar derivative self- interactions. In this paper, we study 
the connection between self-acceleration and the presence of ghosts for a quite generic class 
of theories that modify gravity in the infrared. These theories are defined as those that 
at distances shorter than cosmological, reduce to a certain generalization of the DGP 4D 
effective theory. We argue that for infrared modifications of GR locally due to a universally 
coupled scalar, our generalization is the only one that allows for a robust implementation 
of the Vainshtein effect — the decoupling of the scalar from matter in gravitationally bound 
systems — necessary to recover agreement with solar system tests. Our generalization involves 
an internal "galilean" invariance, under which 7r's gradient shifts by a constant. This sym- 
metry constrains the structure of the tt Lagrangian so much so that in 4D there exist only 
five terms that can yield sizable non-linearities without introducing ghosts. We show that 
for such theories in fact there are "self-accelerating" deSitter solutions with no ghost-like 
instabilities. In the presence of compact sources, these solutions can support spherically 
symmetric, Vainshtein-like non-linear perturbations that are also stable against small fluc- 
tuations. We investigate a possible infrared completion of these theories at scales of order of 
the Hubble horizon, and larger. There are however some features of our theories, that may 
constitute a problem at the theoretical or phenomenological level: the presence of superlumi- 
nal excitations; the extreme sMfr-luminality of other excitations, which makes the quasi-static 
approximation for certain solar-system observables unreliable due to Cherenkov emission; 
the very low strong- interaction scale for txtx scatterings. 



Introduction 



The beauty of General Relativity (GR) lies in the simplicity of its principles and in its empirical 
adequacy in a broad range of phenomena. The modern viewpoint is that GR is an effective 
theory valid below some ultraviolet cut-off, presumably of order Mp ~ 10 19 GeV, and that its 
validity becomes increasingly better the lower the energy, that is the longer the distance scales. 
In particular this viewpoint explains the adequacy of GR at macroscopic and solar system scales 
and implies also its perfect adequacy to describe cosmology at late enough times (and at present). 
However the undeniable beauty of GR has been stained, from almost the beginning, by the mystery 
of the cosmological constant. The inability to provide, to date, a plausible field theory mechanism 
for the absence or smallness of the cosmological constant, perhaps indicates that something is 
missing in the simple infrared picture of GR. Indeed, the strong evidence accumulated in recent 
years that the universe is undergoing a phase of accelerated expansion [1] provides extra motivation 
to consider a modification of gravitational dynamics on cosmological scales. Could late time 
acceleration follow from the very dynamical mechanism that screens the cosmological constant? 
At the moment we have no concrete example in that sense. On the other hand there is instead 
a plausible, and even remarkable, picture where GR is the right infrared decription and where 
the size of the cosmological constant is just explained by anthropic selection [2J. Nonetheless 
the big effort of recent years in trying to come up with, and to understand, modifications of 
gravity has produced an interesting know-how. We have now several examples of theories that 
modify gravity [3j HJ EE E]. In a broad technical sense we may define a modification of gravity 
as a field theory possessing solutions over which new degrees of freedom affect the propagation 
of gravity while the background is not producing any sizable energy momentum tensor [Tj. In 
other words in these theories the graviton propagates differently even in flat space. Examples 
of this type based on purely 4-dimensional field theory are the theories of massive gravity, like 
Fierz-Pauli or more potentially interesting models based on the mechanism of ghost condensation 
[3]. A perhaps more remarkable example is however represented by the Dvali-Gabadadze-Porrati 
(DGP) model [6], based on 5-dimensional spacetime. The DGP model is interesting for several 
reasons. It cannot truly be assimilated to a fluid, as the infrared dynamics involves a continuum 
of states, corresponding to space-time being 5-dimensional, not just a few phonons. It potentially 
represents a stepping stone towards the realization of the old idea of diluting the cosmological 
constant into extra dimensions [HJ, [9] . It provides an alternative explanation for the observed late 
time acceleration of the universe [ID] . Indeed it nicely illustrates that in order to curve space- 
time at long distances it is not necessary to add a zero derivative term, that is the cosmological 
constant, while it is 'enough' a one-derivative term. The latter is in DGP formally represented 
by the 5 dimensional kinetic term, which is more 'relevant' in the tecnical sense than the 4-D 
Einstein term localized at the boundary. Yet another reason of interest in the DGP model lies 
in its successful implementation of the Vainshtein mechanism [11] for decoupling the additional 
modes from gravitational dynamics at sub-cosmological distances. The problem at hand, which 
was first pointed out by van Dam, Veltman and Zacharov [12], is tied to the fact that there is no 
such thing as a purely IR modification of gravity: degrees of freedom with mass of order of the 
Hubble parameter, can expectedly play a role at shorter distances. Indeed, in both Fierz-Pauli 
and DGP these short distance effects are accurately described just by one additional scalar degree 
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of freedom [131 [H], normally indicated by n. The 7r modifies at 0(1) the linearized gravitational 
potential produced by any energy momentum source. Tests of GR within the solar system rule out 
such 0(1) effects, and seemingly the models that produce them. As first remarked by Vainshtein 
long ago, non-linear effects in the dynamics of it are however important. Now, it turns out that 
in the case of Fierz Pauli this does not lead to any benefits, as the non linearities always imply 
the presence of a light ghost state and the consequent breakdown of effective field theory [T5| [T6j. 
However in the case of DGP, under an intersting range of conditions on the energy momentum 
tensor, the Vainshtein effect screens the contribution of 7r below the experimentally acceptable level 
[T71 fT8l fl9] without introducing extra ghost states. Unfortunately the interesting self-accelerating 
solution is not in this class, in that the n is a ghost already on unperturbed de Sitter [T4J 1201 |2T]. 
Now, one crucial property of the 7r Lagrangian C„ in DGP, the one that gives the Vainshtein 
effect a chance to work, is that, while involving higher derivative interactions, it nonetheless gives 
equations of motion of second order in the derivatives [20]. Indeed all of DGP physics at distances 
shorter than H' 1 , including the existence of the self-accelerating solution, is correctly reproduced 
by just the n field coupled to 4D gravity, without any reference to its 5D origin. These facts 
motivate a classification and study of all the purely 4D lagrangians involving just one scalar field 
7r coupled to gravity and to itself in the same peculiar way as in DGP. In particular this approach 
will be enough to study the possible correlation between self-acceleration and existence of ghosts, 
or else to find self-accelerating solutions that lack ghosts. This will be the main goal of this paper. 



1 Local analysis 

1.1 Cosmology as a weak gravitational field 

A FRW universe is a highly non-linear solution of Einstein's equations. Still, at short scales — 
shorter than the Hubble radius — a cosmological solution can certainly be thought of as a small 
perturbation of Minkowski space. For instance, given a comoving observer sitting at the origin, 
we can choose coordinates that for all times make any deviation from Minkowski metric weighted 
by ~ H 2 \x\ 2 , where \x\ is the physical distance away from the observer. Such a coordinate choice 
is not unique. A particularly convenient one is Newtonian gauge, for which g 0i vanishes and the 
spatial metric is isotropic, 

ds 2 = -(1 + 2$)rft 2 + (1 -2^)dx 2 . (1) 
More explicitly, if r and y are FRW comoving coordinates, 

ds 2 = -dr 2 + a 2 (r)dy 2 , (2) 

we perform the coordinate change 

r = t-\H\x\\ y = ~[l + ±# 2 |f| 2 ] , (3) 

where a and H are evaluated at t, rather than at r. Then the metric locally (i.e., for \x\ <C H^ 1 ) 
becomes 

ds 2 ~ - [l - (ff + H 2 ) \x\ 2 ]dt 2 + [1 - \H 2 \x \ 2 ]dx 2 , (4) 
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which corresponds to Newtonian gauge with potentials 

$ = -l[H + H 2 )\x\ 2 , ty = \H 2 \x\ 2 . (5) 

Corrections to the above metric are of order H A x 4 . Notice that Newtonian gauge is a complete 
gauge-fixing only for k ^ (see e.g. ref. [21]); on the other hand our potentials are supported at 
k = exclusively, \x\ 2 ~ 5^(k). We will come back to the resulting extra gauge-freedom later. 
The above form of the metric, eq. (jl]) is convenient because of its being locally a small deviation 
from Minkowski space for all times. That is all along the worldline of the observer sitting at x = 0, 
the metric is 77^ and the connection vanishes. Therefore our gauge-fixing corresponds to Fermi 
normal coordinates. 

Despite the name, Newtonian gauge does not imply any Newtonian limit, and is still fully 
relativistic. However, for a matter dominated universe, all the non-trivial cosmological dynamics 
is captured in our gauge (HI) by Newtonian local physics. Indeed, the Newtonian potential $ above 
obeys Poisson's equation 

V 2 $ = 4vrGp , (6) 

provided that H 2 and H satisfy Friedmann equations. The local fluid velocity field is given by 
the Hubble flow, 

v = Hx, (7) 

which solves Newton's law, 

tf=-V$ = -f. (8) 
a 

The covariant conservation of the stress-energy tensor follows from the Newtonian continuity 
equation, 

p + V ■ (pv) = p + 3Hp = . (9) 

In fact, although we deduced the local form of a cosmological solution from its exact, global 
FRW geometry, for a matter-dominated universe one can derive the correct FRW dynamics just 
from local Newtonian physics. Indeed, assuming isotropy around the origin, Poisson's equation 
with constant p yields $ = 2j ^-p\x\ 2 . The Hubble flow v oc x is the only velocity field that is 
homogeneous — it is invariant under the combined action of a translation and a (time-dependent) 
Galilean boost 0. Newton's law (or Euler's equation) v = — V$ relates the proportionality factor 
(i.e., the Hubble parameter H(t)) to p, via a differential equation. The continuity equation (JHJ) 
then gives the correct time-dependence of p. The fact that this Newtonian derivation gives the 
correct GR results is no accident: for non-relativistic sources and at distances much shorter than 
the curvature radius, the dynamics must reduce to Newtonian gravitational physics. On the other 
hand for a relativistic fluid with sizable pressure, the local interpretation is subtler because of 
relativistic effects, and the coupling to gravity is not fully captured by the Newtonian potential. 

1 Notice that, by performing a time-dependent boost we are in fact going to an accelerated reference frame. 
Since in Newtonian physics acceleration is physical, strictly speaking the velocity field is not homogenous — the only 
inertial comoving observer is that sitting at the origin. However, even in Newtonian physics the only observable 
accelerations are the relative ones — that is the essence of the (Newtonian) equivalence principle. From this viewpoint 
all our Newtonian comoving observers are equivalent. 
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For our purposes it will be more convenient to employ the residual gauge- freedom of Newtonian 
gauge at k = and work with slightly different coordinates, where the metric is locally a small 
conformal deformation of Minkowski space. The downside is that in such a form the metric is 
not near-Minkowskian for all times, i.e. deviations from Minkowski space now will only be small 
locally in space and time. Around x = 0,t = the needed infinitesimal gauge transformation is 

£o = -\ (2H + H 2 ) {\x\H + ft 3 ) , & = \ (2H + H 2 ) Xl t 2 , (10) 

under which the metric (HI) transforms to 



ds 2 



\-lH 2 \x\ 2 + \(2H + H 2 )t 2 (-dt 2 + dx 2 ). (11) 



At the order we are working at H 2 and H can be evaluated at t = and treated as constants. 
The Newtonian potentials now are 

$ = -\H 2 \x\ 2 + \{2H + H 2 )t 2 , ^ = -$. (12) 
1.2 Mixing with it 

Let us consider a long-distance modification of GR where, at least locally and for weak gravita- 
tional fields, deviations from GR are due to a light scalar degree of freedom 7r kinetically mixed 
with the metric. At quadratic order, this amounts to replacing the Einstein action yj—g R with 
sf—g (1 — 2tt)R plus whatever dynamics it may have on its own. For reasons that will become 
clear in the following, even for weak gravitational fields and small ir values we allow for generic 
non-linearities in the 7r sector, whereas we treat the gravitational field and 7r's contribution to it 
at linear order (in the e.o.m.). In this conformal frame we do not couple 7r directly to matter, 
which we assume to be minimally coupled to the metric g^. This ensures that introducing 7r does 
not impair the universality of gravitational interactions. If we demix ir and g^ v by performing the 
Weyl rescaling 

= + 2tc r]^ , (13) 

the action becomes 

S 



J d^x [±M 2 iy ^lR + \ VT"" + C n + ttT%] , (14) 



where once again we are working at quadratic order in h^, and C n encodes the so far unspecified 
dynamics of tc. Now, according to the discussion in the Introduction, if we want the modification 
of gravity mediated by 7r to be an interesting one, 7r's effect on the geometry perceived by matter 
should be dominated by the coupling nT^^ rather than by 7r's contribution to the stress-energy 
tensor sourcing h^ v . For this reason, in this conformal frame we neglect the gravitational back- 
reaction associated with tc (we will come back to the validity of this approximation in sect. 13. 2p . 
In this case, the field equations for h^ v are the same as in pure GR, and so for given sources and 
boundary conditions, their solution will be the GR one. This means that if we know the 'physical' 
metric for given T^ v in our modified theory, we can get the 7r configuration associated with 
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it by subtracting the corresponding GR solution h^ u with the same T^, via eq ([TBI) . However in 
comparing the two metrics we must allow for gauge transformations, because the Weyl rescaling 
eq ([TBI) does not preserve a generic gauge-fixing. The first step is to fix gauges for h^ u and 
such that their difference is conformally flat (if this is not possible, it means that the modification 
of gravity under study is not due to a scalar kinetically mixed with the metric). Then 

V^if "^{J^PlV h/xii ) conformal • (^^) 

However such a gauge choice is not unique. Namely, we are still free to perform infinitesimal 
conformal transformations, under which the metric changes by a conformal factor. Infinitesimal 
dilations shift n by a constant, whereas infinitesimal special conformal transformations add a 
linear piece to it, 

7T -> 7T + C + . (16) 

This means that without further knowledge of the 7r dynamics, 7r is unambiguously determined by 
the geometry only up to a constant and a linear piece. Alternatively, we could try to determine 7r 
by comparing the curvature tensors, which characterize the geometry in a gauge-invariant fashion. 
However at linear order we have 

Rpv = Rpu - 2d li d v ir - rj^Dn , (17) 

which is only sensitive to the second derivatives of it. Once again, the constant and linear pieces 
in 7T are left undetermined. 

Consider now a cosmological solution in such a modified-gravity theory. According to our 
discussion of sect. II. 1[ given the solution a(r) we can choose coordinates such that locally around 
x = 0, t — the metric is a small conformal perturbation of Minkowski space, eq. (fTTj) . If we were 
to do same in GR, at the same value of p and p we would get a different conformal factor, with H 
and H in place of H and H. Then, to extract the 7r configuration that is locally responsible for 
the modified cosmological behavior, we just have to use eqs. (fTTl [T5l) . In terms of the Newtonian 
potential in conformal gauge 

n(x,t) = $-$ (18) 
= -\{H 2 - H%)\x\* + l[2(H-H ) + (H 2 -H*)]t 2 + c + bt, (19) 

where the GR quantities H and H are to be computed at the same values of p and p as H and 
H. We have included the undetermined constant and linear pieces discussed above. In particular, 
given isotropy, the linear piece can only be proportional to t. 



1.3 An example: DGP cosmology 

As an example of our scenario and of the effectiveness of our method, we consider cosmological 
solutions in the DGP model. We will see that locally the modifications to the cosmological 
dynamics are fully captured by a 4D scalar. The four-dimensional modified Friedmann equation 
reads [TU] 

H 2 = ^p + emH = H$ + emH , (20) 
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where Hq = -^rp is the Hubble rate we would have in GR at the same value of p and p, 
m = 1/Ldgp is the critical scale, and e is the sign that distinguishes between the two branches of 
solutions. We can absorb e into the sign of m. Then, the self-accelerated branch corresponds to 
positive m. Taking the time derivative of the above equation, and using the covariant conservation 
of energy, p = —3H(p + p), we get the second Friedmann equation, 

H = -2MfS p + Py( 1 - m l 2H Y 1 = Ho + \m{H/H) , (21) 

where once again H = — ^mt^P +P) i s the corresponding GR value for the given sources. Now 
we want to interpret the deviations from GR as due, at least locally in space and time, to a scalar 
field 7r kinetically mixed with the metric. Then eq. (fl9|) tells us the solution for tt, 

vr(f , t) = -\mH ■ \x\ 2 + \m{H/H + H) ■ t 2 + c + bt . (22) 

We claim that this is the brane-bending mode identified in ref. |14j . To see this, consider the 
e.o.m. that comes from varying the 4D effective action of ref. [2] with respect to tt, 

6Dtt - -1(0 0^)2 + -L(D7r) 2 = -T% = p - 3p . (24) 

If we plug the it configuration (|22|) into this e.o.m., we get a linear combination of the two modified 
Friedmann equations above, namely 12 x (1201) +6 x (121 p . This shows that, despite being intrinsically 
a five- dimensional model, DGP is well described at short distances by a local 4D theory of a scalar 
kinetically mixed with the metric [TU [20]. This description is locally exact, with no need of subtle 
'decoupling' limits — it yields the correct modification to the Friedmann equation as implied by 
the full 5D theory. 



1.4 The structure of C n 

We should now characterize the general structure of C n that corresponds to a potentially realistic 
modification of gravity. On one side we want tt to give rise to an 0(1) modification of the Hubble 
flow, meaning that tt exchange has gravitational strength on cosmological scales. On the other we 
cannot tolerate deviations larger than O(10 -3 ) from Einsteinian gravity at solar system distances. 
It is obviously impossible to satisfy these requests when the tt dynamics is linear throughout the 
universe. However if non-linearities are important, then the extrapolation between different length 
scales is non-trivial and it would be in principle possible to have a viable theory. The DGP model 
with non linearities corresponding to the Vainshtein effect represents an interesting example in 
that sense. So a first necessary request we want to make is that, even in the assumed regime 

2 The present tt has a different normalization from that of refs. [HI [20]. For instance, with respect to ref. [20] 

= 2 • 7r hcre , Mf ere = ±M£r • (23) 
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where tt can be treated as a small perturbation of the geometry, its dynamics is nonetheless non- 
linear. Consider then the general cosmological solution for tt. At subhorizon distances the leading 
behaviour is captured by a quadratic approximation 

tt c = C + B^ + A^xPa? + O (x 3 H 3 ) . (25) 

We remark that the quadratic approximation to the solution is invariant under the combined 
action of a spacetime translation and a shift tt(x) — > tt(x) + b^x^ + c. More precisely 

tt c (x) = tt c (x + 5) - 2A^x v - B^ , (26) 

that is 6 M = — 2A fJlv 5 v and c = —B^. Barring accidents, a symmetry in a solution follows from 
an invariance of the equations of motion and thus constrains the form of the Lagrangian. In the 
case at hand, taking for granted translational invariance, we would deduce that the lagrangian we 
want must be invariant under the shift [f| 

7r(x) -»• 7r(x) + bpX* + c . (27) 

However before taking this deductive step, we have to be more careful, as our symmetry just 
trivially follows from our quadratic approximation: the solution to any scalar field theory can be 
approximated by its Taylor expansion at second order in the coordinates around one point! In 
fact, normally, approximating the solution to quadratic order in the coordinates corresponds to 
linearizing the equations of motion in the field. This is easily seen for an ordinary two derivative 
lagrangian — \{dTr) 2 — V(tt) where the equation of motion linearized in (p = ir(x) — 7r(0) is 

Utp = V'{ir(0)) = const . (28) 

Notice indeed that the above equation admits eq. f[2"5j) among its solutions. Moreover eq. (J2SJ) is, 
as it should be expected, invariant under the shift <p(x) — > (f(x) + b^ + c. Therefore in ordinary 
cases the symmetry n(x) — > ir{x) + b^x^ + c of the approximated solution corresponds trivially 
to an accidental symmetry of the linearized equations of motion. However, as explained above, 
we are here interested in the cases where the tt dynamics, even at local quadratic order in the 
coordinates, is described by a non-linear equation of motion. As follows from our discussion we 
are thus lead to make a second non trivial assumption, namely that the equations of motion of 7r 
be invariant under the shift symmetry fT27j) . This is equivalent to demanding that each tt in the 
equation of motion be acted upon by at least two derivatives. In fact in order to avoid trouble 
with ghost degrees of freedom, a third final request we should make is that the tt equation of 
motion be just of second order. Again the need for this request is exemplified by comparing the 
DGP model to Fierz-Pauli massive gravity. In the latter, unlike the first, the analogue of the tt 
dynamics leads to 4th order equations [15], with a resulting ghost at the onset of non-linearity. 
Our requests imply the equation of motion to take the form 

^ = FiPM (29) 

OTT 

3 Strictly speaking, for a given a solution, characterized by A^ v and B^, combining the transformation in eq. (126[) 
with a translation we obtain a 4-parameter family of transformations where and c are functions of a 4-vector S. 
However the commutator of these tranformations with translations is a constant shift tt — > n + c, so we recover the 
full 5 parameter group. 
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where F is a non-linear Lorentz invariant function of the tensor d^d^n. In the next paragraphs we 
shall classify the Lorentz invariants satisfying eq. fl29|) . It turns out that there is only a handful of 
invariants (more precisely 5 of them in 4D) that do so. Such very constrained class of Lagrangians 
represents a generalization of the ir Lagrangian describing modifications of gravity at sub-Hubble 
distances in the DGP model [Hj. This is indeed not surprising, since our three constraints 
correspond to the three main properties that guarantee viability of the DGP model at solar 
system distances. To conclude, the shift corresponding to c makes ir a Goldstone boson. However 
the vectorial parameter corresponds to constant shift of the gradient <9 M 7r — > d^ir + b^. That is 
the space-time generalization of the Galilean symmetry x — > x + v of non-relativistic mechanics 
(0+1 field theory). We are thus naturally led to term n — > 7r + bx Galilean transformation and to 
call Ti galileon. 

Before proceeding we would like to mention that Galilean invariance of C n can be motivated 
in an alternative and perhaps more direct way. Indeed as discussed in section [L~2l the modification 
brought by ir to the geometry is invariant under the change in eq. ffToT) . That is to say that the 
correction to the Ricci (and Riemann) tensor due to n, AR^ U = —2d^d u 7i — r]^ u Dii, is unaffected 
by eq. (|T6l) . Now, a potentially realistic modification of gravity can be defined as one that corrects 
curvature at 0(H 2 ) at cosmological scales but which is suppressed with respect to ordinary GR 
at shorter scales where the curvature due to local sources is bigger than H 2 . As we said, such a 
behaviour may in principle arise provided the 7r dynamics is non-linear. Then if we further assume 
the ir equations of motion only depend on AR^ V 



we have that corrections to the geometry are directly determined by the distribution of mass in 
the universe (and by boundary conditions). In this situation, we have a better chance to have a 
robust implementation of the Vainshtein effect where AR^ is screened sufficiently near localized 
sources. Thus one viewpoint on Galilean invariance is that it is needed to have a robust and 
economic implementation of the Vainshtein effect. 

We are thus led to classify the Lorentz invariant operators satisfying eq. f[2"§l) . These Lagrangian 
terms must necessarily involve n powers of it and 2n — 2 spacetime derivatives. Of course the 
generic term of this form is not Galilean invariant and correspondingly its contribution to the 
e.o.m. involves also first derivatives of ir. It thus requires some work to find these peculiar invariants 
at each order in n. The simplest possibility is given by n = 1, where we have the 'tadpole' Ci = ir. 
Notice that the tadpole is invariant under Galilean transformations only up to a total derivative 
SCi = b p x p + c = (9 M (6 p x p x M /5 + cx M /4). But this ensures anyway invariance of the equations of 
motion and conservation of a local current. The next to simplest example is the standard kinetic 
energy, (dn) 2 , whose variation under a Galileo boost is 



which is a total derivative. A less trivial example is the DGP cubic interaction of ref. [T4] , \3ir (chr) 2 , 
whose variation is 




(30) 




(31) 



5(Dn(dn) 2 ) = -Ac^d^nDn = -Ad 1 d a [d a ird^ir 



^(dn) 2 ] . 



(32) 
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Again, as follows from Galilean invariance, the variation of the above term w.r.t. tt only depends 
on second derivatives of tt. We can go on and study higher order invariants. In the Appendix 
we show that formally there exists one and only one such Galilean-invariant at each order in 
tt. However, since each derivative Lagrangian term we will construct will be associated with 
one Cayley invariant of the matrix d^d^TT, we will have as many non-trivial derivative Galilean 
invariants as the rank of d^dyTT, that is the number of spacetime dimensions. In particular, in 4D 
we only have five invariants: from the tadpole term, to a quintic derivative interaction @. 

At n-th order in tt, the generic structure of a Galileo-invariant Lagrangian term is (d 2 Tr) n ~ 2 d7rdTT. 
For notational convenience, let's denote by II the matrix of second derivatives of tt, IP„ = d^d v 7T. 
Also the brackets [. . .] stand for the trace operator, and the '•' stands for the standard Lorentz- 
invariant contraction of indices. So, for instance 

[II] Ott ■ Ott = Dtt d^TT . (33) 

Then, up to fifth order in tt the Galileo-invariant terms are 

C x = tt (34) 
C 2 = -I Ott -Ott (35) 
£ 3 = -i[II]07r -Ott (36) 

£ 4 = -l([U} 2 dTT-dTT-2[U}dTT-U-dTT-[U 2 }dTT-dTT + 2dTT-U 2 -dTT) (37) 

£ 5 = -|([n] 3 Ott -Ott- 3[n] 2 8tt ■ II ■ 8tt - 3[n][II 2 ] Ott ■ 8tt + 6 [II] Ott ■ IT 2 ■ Ott 

+2[II 3 ] 8tt ■ 8tt + 3[II 2 ] dn ■ U ■ 8tt - 6 Ott ■ IT 3 • dn) (38) 

The overall normalizations have been chosen to have simple normalizations in the equations of 
motion, see below. Higher order Galileo-invariants are trivial in 4D, being just total derivatives. 

For our purposes it is more convenient to work directly at the level of the equations of motion, 
where there are no integration-by-parts ambiguities, and only second derivatives appear. Defining 
Si = we get 

E x = 1 (39) 
£ 2 = Ott (40) 
S 3 = (Dtt) 2 - (O^tt) 2 (41) 
S A = (□7r) 3 -3D7r(^7r) 2 + 2(^7r) 3 (42) 
£ 5 = (□7r) 4 -6(D7r) 2 (^7r) 2 ^^ (43) 

where by (d fJi d v Tr) n we denote the cyclic contraction, {d^d u Tr) n = [l~l n ]. 

The complete Lagrangian for tt is a linear combination of the above invariants 



C ir = y %2a£ i , (44) 



4 In ID, i.e. mechanics, we just have two: the Galilean kinetic energy ^mi 2 and the linear potential x. 
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where the Cj's are generic coefficients. Likewise, the equation of motion is 



^^ = £ c ^ = - r v ( 45 ) 
i=i 

1.5 Technical analysis 

We now want to investigate the existence of well-behaved self-accelerating solutions. By 'self- 
accelerating' we mean solutions for it that correspond to a deSitter geometry even for vanishing 
stress-energy tensor. Eq. (fl9|) gives us the relevant tt configuration, 

TTdsOr) = -\H 2 x^ , fy&TTds = -\H\ U , (46) 

where we used H = H = 0, because the stress-energy tensor vanishes, and H = 0, because we are 
looking for a deSitter solution. The curvature H 2 is determined by the equation of motion ff4"5]) . 
with T 7 ^ = 0. Indeed plugging the deSitter ansatz into the e.o.m. we get an algebraic equation 
for H 2 , 

ci - 2c 2 H 2 + 3c 3 H 4 - 3c 4 H 6 + |c 5 if 8 = . (47) 

Then, assuming we have a deSitter solution, we are interested in studying its stability and in 
general the dynamics of local perturbations that do not change the asymptotics. With an abuse 
of notation, we will refer to these perturbations by tt. That is, we perform the replacement 
- * TTds + 7i everywhere in the Lagrangian and in the e.o.m. The background deSitter solution 
fj46l) is Lorentz-invariant. Then Lorentz symmetry is unbroken, and it must be an invariance of 
the action for the perturbations. Furthermore, the background is invariant under the combined 
action of a spacetime translation and of an internal Galileo-boost, 

-> X" + a" , 7T dS -> 7r dS + \H 2 a^ , (48) 

which must also be a symmetry of the perturbations' dynamics. But in the e.o.m. only second 
derivatives appear, and they are automatically invariant under internal Galileo boosts. Also, 
the background's second derivatives are translationally invariant. This means that the e.o.m. for 
the perturbations is separately invariant under spacetime translations and under internal Galileo- 
boosts, as well as under spacetime Lorentz boosts. That is, apart from boundary terms in the 
action, the spacetime Poincare symmetry and the internal Galileo boosts are unbroken. The 
e.o.m. for the perturbations then must be a linear combination of our Galileo-invariants (T4QH43]), 

5C W 



5n 



= J2 d i S i> ( 49 ) 
dS i=2 



but now with different parameters dj. The latter are linear combinations of the original param- 
eters Ci with independent coefficients. The constant term £ 1 does not appear in fT49|) because 
it is precisely the e.o.m. evaluated on the deSitter background (j4Tj) . which vanishes. Explicit 
expressions for the mapping {q} <-> {di} (i = 2,..., 5) are given in the Appendix. The exact 
relation will not be important for us. Suffice it to say that it is invertible. 
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In conclusion, as long as we are interested in solutions with deSitter asymptotics, we can 
forget about the original Lagrangian coefficients Cj, pretend to be in flat space, and parameterize 
the dynamics of perturbations with arbitrary coefficients <ij. The determination of the deSitter 
Hubble parameter does not impose any constraint on the coefficients df. for a given set d 2 , • • • > ^5 
we can compute the corresponding c 2 , . . . , C5 and by choosing appropriately c\ there always exists 
a solution of (14"T|) with positive H 2 . 

In order for our Galilean model to be a fairly realistic description of our universe, it should 
pass a number of zero-th order tests. It should admit a deSitter solution that is stable against 
small perturbations. Furthermore, about this deSitter configuration there should exist spherically 
symmetric solutions that describe the tt field generated by compact sources. These configurations 
should also be ghost-free. Finally, we will require that small pertubations about these solutions 
do not propogate superluminally. As we will see our model fails this last test. We will briefly 
comment on this in sect. [H 

For spherically symmetric configurations, tt = 7r(r), the e.o.m. drastically simplifies. In partic- 
ular, like in DGP, it becomes an algebraic equation for ir'(r). This is because the tt e.o.m. is the 
divergence of the Noether current associated with the shift symmetry tt — > tt + c. For generic tt 
configurations the current J M involves first as well as second derivatives of tt. Once we take the 
divergence, all third derivatives cancel, thus yielding two-derivative e.o.m. However for a radial 
configuration 7r(r) the current cannot depend on tt"(t) — otherwise the third derivatives could not 
cancel in the e.o.m.: 

W,*V) = = \r + £^(r) + ^"(r) + % . (50) 

Then the current only depends on tt'{t) (and on r), and the e.o.m. is easily integrated via Gauss's 
theorem, thus yielding an algebraic equation for 7r'(r). Moreover, since each Galilean term in 
the e.o.m. is of the form S n+ \ ~ (ddir) n , the corresponding piece in the current is schematically 
Jn+i ~ {ddTr) n ~ l Btt . For radial configurations 8tt ~ tt' and ddir ~ tt", — , and since the current 
cannot contain tt", we have 

Jn+i ~ r ■ (tx'/t) 71 . (51) 
This means that the e.o.m. reduces in fact to an algebraic equation for (jr'/r). Explicitly, we have 

(52) 
(53) 
(54) 
(55) 

The vanishing of the last contribution is due to the fact that for a time-independent tt the matrix 
dydyir has rank three, thus the fourth-order invariant identically vanishes. For a source localized 
at the origin, p = M5 3 (r), 

^ = Y^diSi = ^r 3 [d 2 (TT'/r) + 2d 3 (Tr'/r) 2 + 2d 4 (7r'/r) 3 ] = M5 3 (f) , (56) 



£2 - 


1 

r 2 




£3 - 


2 

r 2 




£4 - 


2 

-> pr 




£5 - 


-> 
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we get 

d 2 (ir'/r) + 2d, (tt'/t) 2 + 2d, (vr'/r) 3 = ^ . (57) 

Existence 

We see that the existence of a radial solution for all r's constrains the coefficients di to some 
extent. Indeed varying r from zero to infinity makes the r.h.s. of eq. (}5"Tj) span all positive real 
numbers. Let's call P(y) the polynomial of y = ir'/r on the l.h.s., 

P{y) = d 2 y + 2d 3 y 2 + 2d 4 y 3 . (58) 

Then, for a solution to exist P(y) must be invertible when it is positive. We already know that 
d 2 must be strictly positive for the background (deSitter) solution to be stable against small 
perturbations: d 2 is the coefficient of — |(<97r) 2 in the perturbations' Lagrangian. Now, far from 
the source y = n'/r is small, hence the linear term in P(y) dominates. Then, at small y's, positive 
P implies positive y. P is invertible if P'(y) does not cross zero. Thus, in moving to larger y's 
(smaller r's) P has to increase. As a consequence, P is positive when y is. This means that the 
requirement that P(y) be invertible whenever it is positive is equivalent to ask that 

P'{y) = d 2 + Ad 3 y + Qd A y 2 >0 My > . (59) 

A radial solution exists iff this condition is obeyed. At large y (small r), the quadratic term in 
P'(y) dominates, thus d± has to be non-negative. Then eq. (|59|) is satisfied for every y > if we 

also have d 3 > — \J\d 2 d 4 . In conclusion, there exist spherically simmetric solutions with deSitter 
asymptotics iff 

[ d 2 >0 

(60) 




| d 2 G?4. 



Stability 

We now move to study the stability of the radial solution 7r (r) against small perturbations, 7r = 
7r (r) +<f. Given the symmetries of the solution, the quadratic Lagrangian for small perturbations 
will have the form 

S v = \j d A x [K t {r){d tV f - K r (r)(d r cf) 2 - K n (r)(d n ^) 2 ] , (61) 

where (dci(p) 2 is the angular part of (Vy?) 2 , and the kinetic coefficients Ki depend on r through 
the background radial solution 7To(r) and its derivatives. That fact that the full e.o.m., eq. fH9|) . 
only contains second derivatives guarantees that small perturbations about any solution will have 
a two-derivative quadratic Lagrangian, despite the formal appearance of many high- derivative 
terms in the original Lagrangian for 7T. For the solution to be stable against short wavelength 
perturbations, all the K^s must be positive (see a related discussion in ref. [20J ) . 
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Rather than obtaining the K^s by expanding the full 7r-Lagrangian at quadratic order in <p, it 
is faster to compute them by expanding the e.o.m. at linear order, 



£[tt + <p] 



5ip 



-K t (r) (p + -d r (r 2 K r (r) d rV ) + K a (r) d&cp 



(62) 



where is the angular part of the Laplacian. In particular, we can get the radial kinetic coefficient 
from the radial e.o.m. Evidently, comparing eq. (I56I) with the radial part of eq. (1621) . 



K r (r)=P'(y) 



y=K/ r 



(63) 



where P(y) is the polynomial defined in eq. fl58l) . The positivity of K r is thus guaranteed by the 
existence of the radial solution itself, eq. fl59|) . 
For the other K^s we get 



K t (r) 
Kn(r) 



1 d 

3r 2 dr 
J_d_ 

2r dr 



r 3 {d 2 + 6d 3 y + 18d 4 y 2 + 24d 5 y 3 ) y=</r 
■ 2 {d 2 + Ad 3 y + Qd iy 2 ) y=< ^ 



(64) 
(65) 



These depend on y as well as on its derivative with respect to r. Also they explicitly involve r. 
Thus it appears that the question of their positivity is now a differential problem rather than an 
algebraic one. However we can eliminate y'(r) in favor of y and r through the implicit function 
theorem: we know that the background solution satifies 



Then, on the solution 



F(r,y) = P(y) 



y'(r) 



M 
Airr 3 



d r F 



. 



d y F 



This yields 



d\ + 4:d 2 d 3 y + 12(d§ - d 2 d 4 ) y 2 + 24(rf 3 rf 4 - 2d 2 d 5 ) y 3 + 12(3d 2 - 4d 3 d 5 ) y 4 



d 2 + Ad 3 y + 6g? 4 y 1 



Kc 



d\ + 2d 2 d 3 y + (Adj - 6d 2 d A ) y 2 
d 2 + Ad 3 y + 6d 4 y 2 



(66) 
(67) 

(68) 
(69) 



Notice that the explicit r-dependence canceled out. This is because both the background e.o.m. 
and the kinetic coefficients K t for the fluctuation only involve the matrix of second derivatives of 
7r. For radial configurations we have ddn ~ it", — . Moreover the background e.o.m. are linear in 
tt", being the divergence of a current only involving 7r'. As a consequence on the solution it" is a 
rational function of tt' /r, which makes the kinetic coefficients Ki rational functions of tt' /r = y. 

Recall that the solution spans all positive values of y. Thus, we have to require that K t and Kn 
are positive for all positive y's. Given the constraints (|60|) . the denominators are automatically 

positive. The „ U _ in K n the, is pos.tive for * > J^ t which mak es the nu_ in 
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K t also positive provided that cfe < In summary, the radial solution with deSitter asymptotics 
exists and is stable against small perturbations iff 



r d 2 > 

d 4 > 

4 > v 

< ^ 

"5 ^ 4 <f 3 



| 0?2 <^4 



(70) 



Increasing further our level of pickiness, we can ask whether small fluctuations of our solution 
propagate at subluminal speeds. The speed of fluctuations propagating radially is simply 



K r 
K~t 



(71) 



Given the constraints imposed by existence and stability, eq. (I7D1 . this is forced to be superluminal 
at large distances from the source (small y). Indeed 



l + A^y + 0(y 2 )>l. 



(72) 



Notice that the sign of the deviation from the speed of light is positive because of the requirement 
of existence and stability of spherical solutions, independently from the presence or not of the 
self-accelerating background. In moving to smaller distances (larger y) the propagation speed can 
be made subluminal. This can be easily seen from the difference 

K r -K t = Ad 2 d 3 y + 4(4 + 6d 2 d A ) y 2 + 2i(d 3 d 4 + 2d 2 d 5 ) y 3 + A8d 3 d 5 y 4 , (73) 

which becomes negative at large y if d$ is negative. 

On the other hand, angular excitations are always subluminal: their propagation speed is 

2 _K n 



(74) 



The difference between the numerator and the denominator is, apart from an overall positive 
factor 



K n - K t oc -2d 2 d 3 y - (8dj - 6d 2 d 4 ) y 2 - 2A(d 3 d 4 - 2d 2 d 5 ) y 3 - 12(3^ - M 3 d 5 ) y 4 



(75) 



which is negative, as each power of y has a negative coefficient. We end up with the following 
constraints: 

( d 2 >0 
d 4 >0 

d 3 > xfhhdi (76) 



4 < 
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2 Some phenomenology 



In this section we will broadly discuss the scales associated with classical non-linearities and 
with quantum effects. The analysis can be straighforwardly applied also to the fully conformally 
invariant case we discuss in sect. [31 



2.1 The coupling to matter 

The galileon Lagrangian has the general form 

C = f 2 dirdnF(ddir/H 2 ) (77) 

We have directly chosen the scale controlling higher orders in ddir to be the present Hubble 
parameter H , so that the tt contribution affects at 0(1) the present visible expansion, f 2 is 
an overall coefficient that gives the tt Lagrangian the correct dimensions. We can think of / 
as measuring the coupling between tt quanta and matter. Indeed in terms of the canonically 
normalized tt c , upon demixing at linear order, the coupling to matter is ~ jTT c T^^. The scale / 
is bounded by considering the effects, direct and indirect, of tt on the geometry at scales shorter 
than or comparable to 1/Hq. An upper bound is determined by considering the gravitational 
backreaction to the tt energy momentum tensor. At sub horizon scales we have tt ~ H^x 2 , so that 
eq. (ITT)) gives AT^ U = 0(PHqX 2 ). Requesting that AT^ does not lead to a curvature ^> Hq 
at scales just below Hq 1 implies / < Mpi. A comparable lower bound is instead determined by 
considering the direct contribution of tt to the geometry perceived by matter. A bound is thus 
provided by those gravitationally bound systems where the adequacy of GR has been directly 
tested to work at least at 0(1). To derive a constraint a qualitative estimate of the effect of it 
is sufficient. To do so we must figure out which regime, linear or non-linear, characterizes the 
dynamics of the perturbation it over the cosmo logical background TTds- Very much like in the DGP 
model, and as discussed in section 1.5, tt can be treated as linear whenever ddir < 0(Hq), which is 
normally far enough from a localized gravitational source; while ddir > O(Hq) corresponds to the 
non-linear Vainshtein regime, that normally sets in near enough to sources. To get a qualitative 
idea consider a source of total mass M corresponding to a Newtonian potential ~ j^-M/r = 

Rs/r. At large enough distance where the linear regime applies the galileon contribution to the 
metric goes like tt ~ j^M/r = (Mp/ f) 2 Rs/r so we have 

H 2 ~ f 2 r3 H 2 P H 2 U8J 

where by 1Z we indicate the Riemann curvature induced by the source. Now, a system behaves as 
gravitationally bound when the local gravitational attraction wins over the tidal force associated 
with Hubble expansion, that is when 1Z > Hq. This equation shows that when / < Mp, for 
any gravitationally bound system it is in the Vainshtein regime. To obtain a bound on the case 
/ < Mp we thus have to estimate the effects of it by using the behaviour in the Vainshtein region. 
From equation (jSTj) we see that at short distances the quartic coupling dominates the equation of 
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motion. The coefficients d n according to equation (77) scale like 



dn = fj^h (79) 



with a n ~ 0(1), then the asymptotic behavior in the Vainshtein region is n ~ (MHq/ f 2 ) 1 ^ 3 r. 
The relative size of the 7r-potential and the Newtonian one is 



vr /M P \ 2/3 # 4/ V (M P H^ 2/3 



h N \ f J Rf \ f n 



Although the strongest bound on ft/h^ comes from tests of gravity within the solar system, that 
bound does not translate into the strongest bound on /. This is because 1Z within the solar 
system is much bigger than Hq. The strongest bound on / comes instead from more loosely 
bound systems like galaxy superclusters, at whose scales 1Z is just a bit bigger than Hq, and for 
which general relativity seems valid at 0(1) or better. This translates into a lower bound roughly 

/ > (81) 
a few 

In the end we basically have no choice other than / ~ Mp. The fact that both the upper and 
lower bound on / almost coincide is yet another illustration of how difficult is to come up with 
consistent modifications of gravity. Notice on the other hand that in the DGP model, where 7r 
originates from the 5D geometry, one automatically obtains / = Mp. 



2.2 Strong interaction scales and retardation effects 

With the scale / basically constrained to be Mp we can now discuss the issue of quantum correc- 
tions to the classical picture we developed so far. Our discussion essentially follows section 7 of 
ref. [20], where the same issue was studied in the case of the DGP model. However, as we shall 
see, some important differences arise between the two cases. Using / = Mp in eq. (1771) and going 
to the canonical field n = n/Mp, we have 

C = dTrdifF(ddTf/A 3 ) A 3 = M Pl H 2 , (82) 

corresponding to a 'naive' strong interaction scale A ~ (1000 km) -1 precisely like in DGP. Here 
we say 'naive' because, while this result applies around flat space, we expect modifications around 
non trivial solutions. This is because in the Vainshtein region the propagation of the excitations 
ip = 7r — 7To is dominated by the higher order terms in the lagrangian, thus invalidating the above 
estimate. In order to understand what happens we shall focus on the spherical solution that 
was discussed in section 1.5. In the case of DGP it was shown in ref. [20] that, for a spherical 
solution, inside the Vainstein region if = tt — ttq acquires a large kinetic term. This damps ip's 
quantum fluctuations and raises the effective quantum cut-off to A = A(Rv/r) 3 ^, where the 
Vainshtein radius of a spherical source is Ry = Rs/H^. For instance, on the surface of earth one 
has A^ 1 ~ 1 cm. From the discussion in section [27T1 we have 

Tro-^-r — y=Hl~- — »^ (83) 
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y is the parameter controlling the magnitude of non-linear effects inside the Vainshtein region. For 
instance, from the discussion after eq. (|6T|) . we have that the kinetic coefficients in the radial and 
temporal directions K r and K t generically grow like y 2 . On the other hand generically goes 
like a constant, indicating that angular fluctuations are not generically damped when approaching 
the source. Indeed from eq. (161?1) we see that this constant behaviour arises when 7^ 0, while 
in DGP g?4 = g?5 = and K r , K t , Kq ~ y ^> 1. As we shall see in a moment the lack of 
damping for angular fluctuations is an unpleasant novelty compared to DGP. In order to study 
how the fluctuation <p self-interacts, it is useful to write the effective action by distinguishing x±, 
the spacial coordinates tranverse to Vtio (the angular coordinates for our solution) and x\\ = (t, r). 
Since strong coupling in our derivatively coupled theory is a short distance phenomenon, in order 
to assess what the effective strong scale is at a given point, it is enough to zoom in on that point 
and 'linearize' our coordinates and treat them as cartesian. Moreover, precisely as we did with 
the quadratic lagrangian for <p, we can work out the interaction terms working directly on the 
equations of motion. Qualitatively (that is without paying attention to the different contractions 
of indices) the structure of the effective action is 

S eS = J d\\d 2 x^{d\w) 2 + {d^) 2 + a^{d^) 2 dlv 

+ a^{d^) 2 d 2 ^ + a^{d^fd 2 V + ..] , (84) 

where the dots indicate less important interactions of cubic, and higher, order. Notice that the 
interactions involving only transverse derivatives are not enhanced by y factors. This arises because 
for a linear solution 7r oc r one has didjTT oc 6^ with S^j representing the metric in the transverse 
directions. Then if we focus on terms involving only transverse derivatives we have that both the 
background didjTto and the fluctuation are rank two matrices living in the transverse space. But 
such rank two object vanish beyond cubic order because of the rank condition on our Galilean 
invariants. To estimate the strong scale at a given point, we can treat y as constant. The difficulty 
with the above lagrangian is the lack of isotropy. One way to proceed is to change variables and 
make the kinetic term canonical. This can be achieved by working with rescaled coordinates along 
the transverse direction 

~ d 2 x_\_ 
xx = x±/y — > d± = yd±, d 2 x±_ = — — , (85) 



y 2 



so that the action is rewritten as 



5eff = / G^£||<i 2 X_|_ 



(d^) 2 + (d ± ^) 2 + a 3 ^0 ±V ) 2 d 2 ±V 

+ «4^(^) 2 ^ + a 5 ^(9||^) 2 ^ + ...] . (86) 

In this form we can easily estimate the scattering amplitudes 2^2 that arise working at second 
order in the cubic interaction terms. The most dangerous terms are the first two. The third term 
gives the same cut off as around flat space, that is ~ A. The first cubic interaction, the one 
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involving only transverse derivatives, gives an amplitude that roughly scales like 



y4 p8± 1 A ( o 7) 

2 ^ 2 ~ A6 max(Pl, pjj) y 4 A6 max^i/y 2 , pj) ' { } 

where the 1/p 2 factor is an estimate of the propagator and where we have rewritten the amplitude 
also in terms of the unrescaled momentum along the transverse directions. Notice that we have the 
dispersion relation 00 = pjj = \J {p\) 2 + v\lv 2 ~ P±/V- The above amplitude is clearly maximized 

by scattering waves moving along the transverse directions p\\ ~ p±/y- In that case we deduce an 
effective momentum cut-off 

A P ~^ 3 A= (^) V3 A. (88) 

However because of the anisotropic dispersion relation the corresponding energy cut-off is much 
lower 

A £ ~F 2/3 A= (^j ^ A. (89) 

One similarly estimates the amplitude determined by just the second cubic interation term, the 
one involving both transverse and longitudinal derivatives 

a 2 ^ 2 „ y 2 p4±p * = 1 ^ (90) 

~* A 6 max(p^_,p|) y 2 A 6 max(p 2 /?/ 2 , p|) 

With the above amplitude it is convenient to consider two regions: i) p» ~ p±/y and ii) p» > p±/y. 
In the first region we find the same A p ~ y 2 ^A and A^ ~ y _1 / 3 A, while in the second region we 
find that the cut-off is ~ y 1//3 A for both momentum and energy. Finally the third cubic interaction 
term, the one involving only parallel derivatives, implies a cut-off ~ A in both momentum and 
energy. The combined result of this discussion is therefore that the space cut-off remains the same 
~ A but the energy cut off is significantly lowered due to the scattering of the slow moving modes 
along the transverse directions. To get a rough idea we can for instance apply the above results 
to the earth-moon system, neglecting all other sources of 7r other than the earth. For the earth 
we have Ry ~ 10 19 cm and the radius of the moon orbit around the earth is Rm ~ 4 x 10 10 cm 
and thus on length scales of the order of Rm the energy cut-off in time units is 



2/3 

T E = A- E 1 =A~ 1 (^) ~10 3 s. (91) 



Ry 

Rm 

The domain of applicability of our theory is therefore rather limited. Anyway, if the one we just 
illustrated were the only limitation we would still, at least in principle, be able to defend the 
estimate of the effects of n on the moon orbit, which was one of the interesting features in DGP. 
After all in that case the time scale, a month, is much bigger than T E . However there are worse 
limitations (or better complications) in computing the effects of n within the solar system already 
at the classical level. This is because of the big retardation effects associated with the slow motion 
of the n excitations. In order to compute the 7r field generated by the earth we should first find 
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the 7r solution due to the sun, and then solve the non linear problem around the earth. Now, given 
the Vainshtein radius Rvq associated with the Sun mass and the earth orbit radius Re we have 
that angular fluctuations propagate with a speed 

v „ x c ~ l(T 8 c , (92) 
Rvq 

which is orders of magnitude smaller that the velocity of the earth on its orbit. The computation 
of the 7r field around the earth could then not be performed in the static limit. Correspondingly 
Cerenkov radiation in the transverse direction would instead be emitted due to the motion of the 
earth. The basic conclusion is thus that the lack of enhancement of transverse kinetic terms around 
a spherical source causes trouble. On one side we have a technical trouble at the classical level, as 
huge retardation effects prevent a simple static computation of the background. On the other side 
the have a conceptual trouble at the quantum level, due to a rather low cut off, worse in energy 
than in momentum, again because of retardation. Maybe our conclusion for the solar system are 
too drastic in that we have not considered other effects, like that due to the galaxy, which will 
enhance the transverse kinetic terms above their 0(1) value we used for our study. In principle 
one could try and estimate these effects. For instance by studying the 7r background around two 
separated spherical sources. The lack of spherical symmetry makes this a much tougher problem 
that the spherical one we studied, and is beyond the scope of the present paper. 

2.3 Special cases 

So far we have discussed the most general lagrangian C n in eq. (PHI) , with all the coefficients c, 
different from zero. In the following we examine if restricting to a subclass of these theories can 
alleviate some of the difficulties found in the previous sections. 

From the qualitative analysis of sec. 12.21 we can see that the presence of the higher order 
galilean terms £4,5 is responsible for the slow motion of the 7r excitations that in turn restricts 
the domain of applicability of the theory as one goes to shorter distances. It's interesting then 
to consider what happens if we choose <i 4 = d 5 = 0, or equivalent c 4 = c 5 = 0. In this case the 
only interaction left is the cubic one of the DGP model and the kinetic terms K r , K t and Kq in 
the region where classical non-linearities are important have an isotropic enhancement oc y. This 
large kinetic energy suppresses quantum fluctuations if with the result that the position dependent 
cutoff increases, for example on the surface of the earth we have A^ 1 ~ 1 cm [20] • The difference 
however is that — contrary to DGP — there can still be a self accelerating solution without ghosts 
because of the presence in the theory of the tadpole C\K. Moreover the conditions ([70]) can be 
satisfied when c? 4 = <i 5 = 0. 

The case c\ = could be regarded as theoretically preferable, as long as both the Minkowsky 
and deSitter solutions are stable. For instance in applications to early cosmology this would 

5 Notice in passing that the possibility to choose also dj, — seems to be compatible with the constraints (|70|l. 
What goes wrong in this case is that without any non-linear effects in the dynamics of 7r we can have deSitter 
asymptotics but there are 0(1) deviation from GR already at solar system scales. As we stressed throughout the 
paper, the implementation of the Vainshtein effect is a necessary condition to reconcile theories with a lagrangian 
of the form (|14|) with present astronomical observations. 
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potentially allow to describe inflation and Big Bang cosmology within the same effective theory. 
Moreover for c\ — the deSitter solution would arise in the absence of a potential, which is 
one possible definition of self-acceleration. However one should keep in mind that also the case 
ci 7^ gives a quite novel way of generating a deSitter solution: there is, sure, a potential for the 
scalar ir but, unlike the case of a cosmo logical constant it does not have any, even approximately, 
stationary points! 

Anyway with the above caveat in mind it makes sense to better analize the case C\ = to 
find out if our zeroth order requirements can be satisfied. Can this case be made compatible also 
with the existence of a spherically symmetric solution? Contrary to the analysis carried out in 
sect. 11.51 when c\ = equation ( 14T1) gives a constraint on the coefficients 4: we can rewrite it 
factoring out an overall H 2 and defining Q(^H 2 ) the resulting polynomial: 

Q{y) = 4(4 + 34 y + 64 y 2 + 64 y 3 ) . (93) 

Then, for the background deSitter solution to exist at all, Q(y) must have a zero at positive y, 

3y >0 | Q(y„) = 0. (94) 

If 4 = 4 = we can see immediately that (1941) cannot be satisfied: it's clearly incompatible with 
the simultaneous requirement of stable fluctuations around deSitter (d 2 > 0) and existence of a 
radial solution for all r (da > (|60|) ). This is the well known lesson of DGP. 

However even when 4, 4 ^ there is some tension between the two algebraic eqs. (1591 l94"l) . 
In spite of that, it's easy too see that they can be both satisfied if we add to the conditions ( 1601) 
the requirement 4 < 0. Therefore when c\ = the presence of the quintic term with a negative 
coefficient in the Lagrangian becomes crucial to make the deSitter background and spherically 
symmetric solutions coexist. Notice that 4 < is also compatible with the requirements of 
stability and subluminality of the fluctuations f[T5"j) . at least close to the source. 

We can still take 4 = trying to avoid the problem related to the speed of the angular 
fluctuations. In this situation we have K r , Kq ~ y while K t ~ y 3 . Even before paying attention 
to the UV cut-off, this case looks right away uninteresting since the retardation effects are even 
worse than the general case discussed in the previous section. Indeed all modes have a speed 
v ~ c/y. Now, from eq. (157|) at 4 = 0, at a distance r from a spherical source with Vainshtein 
radius Ry we have y ~ (Ry/r) 3 ^ 2 . So in the case of earth orbiting the sun v ~ 10~ 12 c. 

3 The IR completion 

3.1 Self-acceleration from symmetry 

We now consider the issue of whether there exists a "global completion" of the local 7r dynamics we 
have been discussing so far. It may be that at large distances the description in terms of 7r breaks 
down, like in DGP, where at horizon scales the four-dimensional description is no longer applicable. 
Here we want to entertain the possibility that instead, the physics at large scales is well described 
by a four- dimensional Lagrangian for n. To this purpose, it is instructive to characterize the 
acceleration of the universe in terms of symmetries. If the universe asymptotes to a deSitter phase 
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at late times, it is in fact appoaching an enhanced symmetry configuration. Indeed a generic FRW 
solution is invariant under spatial translations and rotations, whereas deSitter space is maximally 
symmetric, with as many isometries as Minkowski space. The cosmological constant is perhaps 
the most elegant way of realizing such symmetries: having a stress-energy tensor proportional to 
the metric, it does not introduce any new tensor thus being consistent with the highest possible 
number of symmetries a metric may have. It is certainly the most economic way, in that it does 
not introduce any new degree of freedom. A priori, a possible alternative is the usual way in 
which symmetries are realized in field theory: at the level of the Lagrangian, off-shell. That is, if 
a solution to the field equations has some symmetries, barring accidents these must be shared by 
the Lagrangian that gives rise to such a solution. The most economic possibility in this class is to 
realize the desired symmetries introducing the least possible number of new degrees of freedom. 
As we will see shortly, for the deSitter symmetry group one scalar suffices. 

Let's start by considering the limit in which gravity is decoupled Mpi — > oo. Notice that 
the full Poincare-like symmetry group of four- dimensional deSitter space — the generalization of 
translations, rotations, and boosts — is 5*0(4,1), that is, the group of five- dimensional Lorentz 
transformations. Indeed .D- dimensional deSitter space can be embedded into (D + l)-dimensional 
Minkowski space as an hyperboloid, which is a Lorentz-invariant submanifold. This suggests 
a quick way of constructing 4D theories that realize the deSitter symmetry group as an internal 
symmetry: write down the most generic effective Lagrangian for a 4D brane living in 5D Minkowksi 
space. The presence of the brane itself will necessarily break spontaneously some combinations 
of 5D translations, rotations, and boosts. The two largest possible residual groups are the four- 
dimensional deSitter group 50(4, 1) — when the brane is an hyperboloid — and its flat limit, the 
four-dimensional Poincare group. The existence of deSitter solutions will ultimately depend on 
the Lagrangian, but may be generically expected given the symmetry arguments above. From this 
geometrical perspective, it is clear that we are introducing just one new 4D scalar — the position 
of the brane. Notice that in this approach the fifth "extra-dimension" is fictitious — no field 
propagates in it. It is simply a useful geometric picture of the new degree of freedom. However, 
if we take it more seriously, we introduce gravity in the 5D ambient space, and we also gauge 
gravity on the brane, what we get is nothing but the DGP model 0. From this viewpoint it is not 
surprising that DGP admits "self-accelerating" deSitter solutions. 

Another possibility, which is the one we will analyze in the following, is to make the 4D 
Lagrangian conformal invariant. The 4D conformal group is 50(4,2), whose maximal subgroups 
are 50(3,2), four-dimensional Poincare, and 50(4,1) — the isometry groups of four-dimensional 
Anti-deSitter, Minkowski, and deSitter spaces, respectively. It is well known that the spontaneous 
breakdown of the conformal group to any of these subgroups can be achieved by a single scalar, 
the dilaton [23] • Indeed the 4D conformal group is the isometry group of 5D AdS space, and 
the dilaton can be thought of as the position of a 4D brane in 5D AdS. Then each of the above 
breaking patterns corresponds to a maximally symmetric brane geometry. 

Let us stress again that for the moment we are decoupling gravity. We will see in sect. 13.21 
that reintroducing dynamical gravity drastically changes the picture at large distances. However 
for the local analysis we are carrying out, at distances much smaller than the Hubble scale, the 

6 For DGP an orbifold projection is also necessary. 
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effects of dynamical gravity are negligible. Let's therefore consider the case in which matter is 
minimally coupled to the metric 

9^u = e^r]^ , (95) 

where the exponential is just a convenient non-linear completion of eq. (1131) . and we "froze" 
dynamical gravity, g^ v = rj^. Then it non-linearly realizes the conformal group: 

• Dilations 

D : n(x) tt'(x) = tt(Xx) + log A ; (96) 

• Infinitesimal special conformal transformations 

: tt(x) -> 7i'(x) = 7i(x+ [ex 2 - 2{c-x)x)) - 2c^ ; (97) 



• Translations 

: ir(x) — > n'(x) = 7r(x + a) . (98) 

• Boosts 

M„ v : ?r(x) tt'(x) = tt(A • x) . (99) 

The three possible maximally symmetric solutions (AdS, Minkowski, dS) correspond to dif- 
ferent unbroken combinations of the above generators. Clearly, Minkowski space corresponds to 
a trivial (constant) configuration for tt, for which translations and boosts are unbroken, whereas 
dilations and special conformal transformations are spontaneously broken. More interesting are 
the other solutions. The deSitter solution is [25] 



(100) 



l + \H 2 x 2 



which manifestly leaves boosts unbroken. It is also invariant under the combined (infinitesimal) 
action of and P M , 

S^ = P^-\H 2 K^: tt(x) -v ir'(x) = n(x + a- \H 2 (ax 2 -2{a-x)x)) + \H 2 , (101) 

as can be straightforwardly checked. Indeed the four generators together with the boost 
generators form the full algebra of 50(4, 1), the deSitter symmetry group. Likewise, the AdS 
configuration corresponds to 

which preserves boosts as well as the combination = P M + ^k 2 K^. The unbroken symmetries 
make up the AdS symmetry group, SO(3, 2) [25] . 

If we start from a conformally invariant Lagrangian for tt, we may generically expect one of 
these configurations to describe the vacuum. Which one will be a solution to the dynamics, and 
with what curvature scale, clearly depends on the Lagrangian. However, we just have to ask 
whether we can generate all the Galilean-invariant terms we need for a locally deSitter solution to 
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exist, be stable, and admit stable spherically symmetric perturbations, starting from a conformally 
invariant Lagrangian. If we can, then by symmetry the deSitter solution will exist globally, and 
will have everywhere the same local properties (stability, etc.). 

To construct conformally invariant Lagrangians, we simply write down all possible diff-invariant 
Lagrangian terms involving the fictitious metric ([9"5l) . In fact, the conformal group is a subgroup 
of diffeomorphisms — the subgroup that leaves the metric conformally flat. The simplest term is 
obviously 

l 'd l Xy/=g = [ d 4 xe 47T , (103) 



which corresponds, upon linearization in tt to the term proportional to tt in the galileon lagrangian 
of the previous sections. It is clearly conformal invariant, upon changing the integration variable. 
For instance, under special conformal transformations, eq. ( 1971) . we get 

d 4 xe 4 ^ x) -> / ^ e M*+(^ 2 -2(c-z)*)) e -*v* = f d 4 x'e^ {x,) . (104) 



The next terms, in a derivative expansion, are the curvature invariants R, R 2 , R^R^, R tJ, upa R tJi upa ' , 
and so on. Let's look more closely at their structure in terms of tt. The Riemann tensor for the 
conformally flat metric eq. (1951) is schematically 

R? upa ~ dnd-n + d\ . (105) 

Thus a generic n-th order curvature invariant is of the form 

Jd 4 x^R n ~ Jd 4 xe {4 - 2n)n (dndn + d 2 7r) n 

~ j d 4 x e {4 ~ 2n)7r [(d7rd7c) n + ■■■ + OttOtt (d 2 Tr) n - 1 + (<9 2 vr) n ] (106) 

The Galilean invariants we constructed in the previous sections correspond to the second-last term 
in eq. (I106p . The terms with more 7r's and fewer derivatives per tt — those on the left of the Galilean 
one in eq. (1 106[) — are less important at short distances and small values of tt. For instance on the 
deSitter configuration tt ~ —^H 2 x 2 , the relative importance of one such term with respect to the 
Galilean one is 

dndn ( } ' m 

which is negligible at sub-horizon scales. The same conclusion obviously holds for the exponential 
prefactor in eq. (11061) . which at lowest order we can set to one. In fact at short distances and 
at small n- values, conformal invariance reduces to Galilean invariance. This is manifest from 
expanding dilations (eq. (1961) ) and special conformal transformations (eq. (1971) ) at lowest order in 
7r and x. 

However by the same token, at short distances and for small 7r the Galilean term in eq. (11061) 
is negligible with respect to the one on its right, which has one less tt and two derivatives acting 
on every tt. For instance their ratio on the deSitter solution is of order H 2 x 2 . The terms with 
two derivatives per tt are exactly the terms we did not want in sect. 11.41 because they lead to 
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fourth-order equations of motion, and thus to ghosts. Therefore, we need to find combinations 
of curvature invariants for which these unwanted higher-derivative contributions vanish. For such 
combinations, the most important Lagrangian terms at short distances will be the Galilean ones, 
and our local analysis of the previous sections will apply. Since the tuning that sets to zero the 
(d 2 7r) n terms corresponds to removing a degree of freedom, terms of this type (by definition of 
perturbation theory) will consistently be generated at the quantum level only with a coefficient 
with a small enough size that the ghost will be at the cut-off. This issue is discussed in ref. |20j . 
The first curvature invariant is 



J d A x^gR = J d 4 xe 2w [-6(<9tt) 2 - 6 Dtt] . 



'108) 



Notice that here we do not need (and we could not either) cancel the Ok term. Indeed integrating 
by parts 

(109) 



,2tt 



□tt 



2e 2n (dK) 



and so eq. fl 1 8 [) is nothing but the ir kinetic term. 
The next invariant is 



J d 4 x^ [aR 2 + b R^BT + c R^R/** 7 ] . 



(110) 



Here we have a complication. Upon integration by parts, the three terms above yield all the same 
structure in terms of tt 



-gxR 2 ,(R^) 2 ,(R>* ^ 2 



vpa ) 



oc [(<9vr) 4 + 2(<9vr) 2 Ok + (Dtt) 2 ] 



(111) 

where (<97r) 4 = ((dn) 2 ) 2 . Hence, there is no way of tuning the coefficients in eq. ( 11 101) as to make 
the (Ok) 2 term vanish while keeping a non-trivial Galilean term. However this is an accident 
peculiar to 4D. We can ask whether there is a suitable combination of invariants in d dimensions, 
for which the (Ok) 2 vanishes and which is regular and non-trivial in the cZ — ^ 4 limit — a sort of 
derivative with respect to d. In d dimensions the Ricci tensor and scalar read 



R 



R 



(d - 2)d^nd u n - (d - 2)rj^ u (dTx) 2 - (d - 2)8^71 - -q^D 



71 



-2?r 



-e - [(d -l)(d- 2)(<9tt) 2 + 2(d - 1)Dtt] 
(we will not need the Riemann tensor). Then for instance the combination 



(112) 
(113) 



d xsf^g d 



d-4 



RfiuR^ 
d- 1 



R 2 



d d x e {d ~ A)n 



2(d 



(d-1) 2 . 
idn) 2 DK 



(d-2)(3d-4),_ , 2rn (d 



2(d- 1 



-{dKf 



[Uk) 2 



(114) 



is regular for d — > 4, and linearly independent of eq. (11111) . By combining eqs. (II lip and (11141) 
we can tune to zero the (Dtt) 2 while retaining the Galilean interaction (as well as the subleading 
(chr) 4 term). 
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Alternatively, we can shelve for a second the auxiliary metric g^ v = e^rj^, and look for 
invariants in 4D directly in terms of tt. At the four- derivative level (which corresponds to R 2 
etc.), invariance under dilations (eq. fl96l) ) tells us that there is no e nn overall factor. Then, under 
special conformal transformations (eq. (I97p ) the individual terms transform as 

Jd 4 x(d7T) 4 -> Jd 4 X [{dTif - 8c^7r(«97r) 2 ] (115) 

Jd A x{d7T) 2 n7i -> jd 4 x [(dix) 2 nn + Ac^d^^d-n) 2 - Ac^d^Un] (116) 
jd A x{Uii) 2 -> Jd 4 x [(Dtt) 2 + 8c^7rD7r] (117) 

where on the r.h.s.'s we changed integration variable, x M — > s 1 ' + (c^x 2 — 2(c • x)x^). We thus see 
that at the Lagrangian level the only invariant combination is precisely eq. fillip . However, if we 
require that the e.o.m. be invariant, the Lagrangian need only be invariant up to a total derivative. 
Then the (Ott) 2 term is invariant by itself, since d^dn is a total derivative — see e.g. eq. (132]) . 
This means that we can set it to zero while preserving conformal invariance. 

Notice in passing that Lagrangian terms of the form e'- 4-2 ™-' 77 (d 2 ir) m (dir) 2n ~ 2m are automatically 
invariant under dilations, eq. (196]) . whereas as we saw only some very specific combinations of them 
are invariant under special conformal transformations as well. This goes against the common 
expectation that scale invariance alone in enough to guarantee full conformal invariance. In fact, 
the standard arguments do not apply for higher- derivative terms — see e.g. ref. [27J. 

At the next orders, R 3 and i? 4 , there are no complications. Suitable combinations of curvature 
invariants involving just the Ricci tensor and scalar give the Galilean terms we want, while tuning 
the unwanted (Dtt) 3 and (Dvr) 4 terms to zero. At order R 3 we get (after integration by parts) 

V=g [i R 3 + R(R^) 2 - (R^) 3 } = e~ 2 * [24 A + • • • ] , (118) 

where £4 is the fourth-order Galilean invariant, eq. (|37|) . and the dots stand for terms with fewer 
derivative per ir, which are negligible at short distances. Likewise, at order R A 

= e- 47r [i6£ 5 + ...] . (119) 

In conclusion, we can reproduce any Galilean-invariant Lagrangian of the form eq. ( 144]) as 
the short-distance limit of a globally defined, conformally invariant Lagrangian for it. Then, the 
existence of a deSitter solution as well its properties (stability, etc.) can be studied locally, about 
any point, as we did in the previous section. The symmetries of the full Lagrangian guarantee 
that the deSitter solution we find locally can be promoted to a global one, and the isometries of 
deSitter space guarantee that the stability and the other local properties of the solution will be 
the same throughout space. 

3.2 Gravitational back-reaction 

We now consider the (phenomenologically interesting) case where gravity is dynamical, Mpi < 00. 
In our notation, this amounts to changing the metric seen by matter, eq. (195]) into 

9^ = e 2n g^ , (120) 
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R 2 {R^ V ) 2 + J2 R{R^u) 3 + jq{(Riiu) 2 ) 2 - l(R, 
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of which eq. ffl3l) is the small field limit, and making g^ v dynamical. Also we should covariantize 
the tt action, which corresponds to coupling tt to g^ v . As usual there are many inequivalent ways 
to do so, parameterized by non-minimal couplings of tt to curvature invariants. Now, if we want 
our deSitter solution to survive after the coupling to gravity is taken into account, we have to 
require that gravitational interactions do not spoil the conformal invariance of the tt action — so 
that the existence of a deSitter solution is guaranteed by symmetry. But the only way to couple 
a conformal theory to gravity while preserving conformal invariance is imposing that the action 
be invariant under Weyl rescalings of the metric: 

V(^j;(i)=e 2uW v(i (121) 

If we do so, the dilaton of the CFT and the metric have to appear precisely in the combination 
( 11201) . so that a conformal transformation on tt may be reabsorbed into a Weyl transformation of 
g^. This corresponds to gauging the conformal group. In fact, this is how we constructed the 
conformally-invariant action for tt in the first place — by making the conformal group that acts on 
tt a subgroup of diffeomorphisms acting on g^ v . The full action is then the sum of the manifestly 
covariant operators of last section, where now R^ vpa stands for the Riemann tensor of the metric 
f )120p FI As a result, our dilaton tt disappears from the action — nowhere do tt and g^ v appear 
separately from each other. In fact, via a Weyl transformation we can choose 'unitary gauge', 

tt = 0, g^v = g^u- (123) 

The action now is just a higher- derivative action for the metric, coupled to matter. The problem 
with this action, is that the energy scale that "suppresses" the higher-derivative terms is minuscule. 
Indeed, in the presence of an Einstein-Hilbert term M^R, in order to have a deSitter solution 
with curvature ~ Hi, we need either a non-derivative operator 

r2 tt2 



g M^Hq , (124) 

which is not surprising news, or higher-order curvature invariants weighed by inverse powers of 
Hi, 

Mpi R n . (125) 



Hl n 



The latter possibility, which is the non- conventional one, corresponds to a ultraviolet modifica- 
tion of GR at curvatures of order Hi, which, even at the classical level, is clearly ruled out by 
experiments. 



7 The combination in cq. (| 1 14[) is non-singular in d = 4 only for conformally flat metrics, as manifest from its 
l.h.s. This means that the DGP-like term 

2(<97r) 2 D7r+ (<9tt) 4 (122) 

cannot be coupled to gravity in a Weyl invariant fashion. Given that — as we will see shortly — the Weyl invariant 
couplings we are constructing are not viable, we do not investigate this issue further. Also notice that when using 
a generic metric g^v rather than the conformally flat one of last section, we have new independent higher-order 
curvature invariants, involving for instance the Riemann tensor — which for conformally flat metrics is instead 
determined by the Ricci tensor. 
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We are thus led to the conclusion that viable gravitational interactions of ir cannot preserve 
conformal invariance — a conclusion that is also strengthened by Weyl invariance's being generically 
anomalous at the quantum level. As a consequence, the deSitter solution we found for it in general 
will be spoiled when the gravitational backreaction is taken into account. Also notice that the 
problem of having higher- derivative terms for the metric with huge coefficients like in eq. (I125p . 
will be shared by a large class of Lagrangians where 7r has important non-minimal couplings to 
g^ u . The safest possibility from this viewpoint seems to be minimal coupling — it involves the least 
number of derivatives of g^ compatible with diff- invariance. Let's therefore consider the case of 
minimal coupling. This amounts to covariantize the n action just by contracting Lorentz indices 
with g^, replacing ordinary derivatives with covariant ones, and integrating in d A xy/—g. Also, 
we include an Einstein term M^R, to reproduce Newton's law. Scaling the various terms in the 
galilean limit like in eq. (1771) and using the results of section 3.1 for their relation to the other 
terms, the n-th order Lagrangian term is schematically 

£ ~ -^- 2 e^ 2 ^(dndn + ddir) n . (126) 

Now, as we saw, a deSitter configuration for 7r corresponds to eq. (11001) . For cosmology, it is 
more convenient to use coordinates where the solution only depends on (conformal) time, 

*- = --L . (127) 

The relation between this and eq. (UOOp is just a special conformal transformation on 7r combined 
with a translation. However, as we discussed, when gravity is taken into account conformal trans- 
formations are no longer a symmetry. Therefore the two expressions for the deSitter configuration 
are physically inequivalent — i.e. they produce different gravitational fields. If we are looking for 
FRW solutions for the metric, we have to stick to eq. (I127p . The stress-energy tensor scales like 
the Lagrangian (I126p . which for our deSitter configuration gives 
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The curvature perceived by matter is a combination of that due to its conformal coupling to tt, 
and that of g^ u , which in turn is sourced by the 7r stress energy tensor: 

~ dirdn + ddn + R^ (129) 



V 



1 + 



./ 



2 



(130) 



The latter contribution becomes sizable and spoils the deSitter geometry at conformal times 



ws^-s' (131) 

If we could make / smaller and smaller, we could postpone this moment at will. However we have 
already shown that tests of GR at sub cosmo logical scales constrain / to be not significantly smaller 
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than Mp. So we cannot play with / to make the solutions more stable. Anyway, since proper 
time t depends only logarithmically on rj, this means that the geometry will look approximately 
deSitter only for times of order Hubble. That is, gravitational interactions necessarily spoil the 
deSitter solution at Hubble scales. 

4 Discussion and Outlook 

We conclude with a few remarks on our results and on possible future developments. First, let's 
summarize our findings. We have considered theories that modify the infrared behavior of gravity 
where, at least in a local patch smaller than the cosmological horizon: 

(i) The modification of gravitational interactions is due to an effective relativistic scalar degree 
of freedom 7r kinetically mixed with the metric, or equivalently, universally coupled to matter 
via the interaction ttT^^. 

(ii) This scalar degree of freedom generically decouples from matter at shorter scales thanks to 
important derivative self-interactions (Vainshtein effect). |^| 

In this class of theories, we have argued that the dynamics of 7r must enjoy an internal Galilean 
symmetry. In 4D there are only five possible Lagrangian terms with such a symmetry. This 
makes a thorough analysis of very symmetric solutions possible. First, we demanded the existence 
and stability of a deSitter solution in the absence of a cosmological constant — a self-accelerating 
solution a la DGP, but with no ghost instabilities. Then, we imposed that with these deSitter 
asymptotics there exist spherically symmetric Vainshtein-like solutions, describing the tt config- 
uration about compact sources, and that these configurations are also stable against small per- 
turbations. We found some tension, if no contradiction, in satisfying these requirements. It is 
plausible — although we have no proof of this — that demanding the existence and stability of more 
complicated solutions, like e.g. with several compact sources in generic positions, may lead to 
contradictions. However, we also found that small excitations are forced to be superluminal at 
large distances from a compact source, where ir is in the linear regime, even if we don't require 
the existence of self-accelerating solutions. 

We ran into trouble also as a result of the opposite problem: the very slow speed of angular 
fluctuations when approaching the source due to the presence of the galilean terms £4,5. This 
makes a static computation of the background not trustable because of huge retardation effects, 
which are also at the origin of the very low energy cutoff of the theory. To avoid these difficulties 
we considered the subclass of Lagrangians with 04 = 05 = 0. This requirement is still compatible 
with the absence of ghosts around deSitter background and with the existence a stable Vainshtein- 
like solution around compact sources. The kinetic term in the angular directions is now enhanced, 

8 This is the case for instance in Fierz-Pauli massive gravity and in the DGP model, whereas massive gravity 
theories that crucially rely on spontaneous breaking of Lorentz invariance are left out. 

9 Some form of non-linearities are necessary for screening a universally coupled scalar at short scales, thus 
recovering agreement with solar system tests of GR. The only robust alternative to the Vainshtein effect we are 
aware of is the Chameleon mechanism |30j . There, the non-linearities are in the coupling of the scalar to matter, 
rather than in its self-interactions. 
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as a result quantum fluctuations are suppressed and the effective cutoff increases like in DGP 
above the naive estimate A ~ (1000 km) -1 . 

Finally, we considered whether our local galileon tt can be promoted to a global scalar degree 
of freedom, coupled someway to four- dimensional GR. First neglecting gravity, we promoted the 
Poincare group and our local galilean invariance to the conformal group, with tt as the dilaton. 
The deSitter group is one of the maximal subgroups of the conformal group, and the existence of 
deSitter-like solutions for tt may be generically expected by symmetry considerations. We showed 
that our local Galilean invariant terms are in fact the short-distance limit of globally defined, 
conformally-invariant Lagrangian terms for tt. The problem however is in coupling this system 
to gravity. This necessarily breaks conformal invariance, and as a consequence spoils the deSitter 
geometry at times of order of the Hubble time, and longer. Although this is not a disaster — in 
fact this model may well yield a realistic cosmology, for we have no idea of whether our universe's 
future resembles deSitter space — the symmetry motivation for the 'conformal completion' of the 
Galilean symmetry is somewhat lost. 

Let us stress once again that our local analysis of the Galilean model is quite generic — it applies 
to all theories obeying conditions (i) and (ii) above. Therefore, the partial failure of our attempt 
to 'IR-complete' the Galilean model should not be seen as impairing the local analysis. Our local 
analysis applies for instance to the DGP model — as we showed in sect. II .31 DGP at short distances 
corresponds to a Galilean model with non-vanishing C2 and C3 only. There however, it is manifest 
that the IR completion is not a four-dimensional local theory — it is in fact a five- dimensional one, 
and our tt degree of freedom stops making sense at scales of order of the cosmological horizon, 
where its interactions and mixings with the whole continuum of KK modes cannot be neglected 
anymore. We do not know whether this is going to be the case for a generic Galilean model, or 
whether instead there exist purely four- dimensional local IR completions that realize the deSitter 
symmetry group as a subgroup of a larger internal symmetry group, other than the conformal 
one. Perhaps, as we suggest in sect. 13.14 there exists one such theory where the internal symmetry 
group is the 5D Poincare group. This in a sense would yield a purely four- dimensional version of 
DGP, without the extra complications of 5D gravity. 

On the other hand, since DGP locally reduces to a particular Galilean model, perhaps a generic 
one can be IR completed by a generalization of DGP. Possible generalizations include higher-than- 
five dimensional setups[3I], or the standard five- dimensional setup with non-minimal bulk or brane 
operators. In the latter case, we might try adding higher-curvature invariants in the bulk, but in 
order to keep second-order e.o.m. for tt — a crucial feature of our galilean-invariant Lagrangians — 
we are only allowed to use Lovelock combinations. But in 5D we just have one such combination, 
whereas to reproduce a generic Galilean model we need several independent terms. The other 
possibility is to add localized terms on the brane involving powers of the extrinsic curvature. 
Indeed for small gravitational fields the extrinsic curvature is K^ v oc d^d^ir, so such terms will 
directly yield derivative self- interactions for tt. The problem with this is that they involve too 
many derivatives — two per tt field. Another possibility would be to add new five-dimensional 
degrees of freedom besides gravity, or, as we said, to consider higher- dimensional spacetimes. We 
leave this for future investigations. 

We now turn our attention to the issue of the ultraviolet completion. From this viewpoint, 
our model is just a generalization of the 4D effective action of DGP that gets rid of the ghost 
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on the self-accelerated branch, but that nevertheless shares with DGP the other pathologies iden- 
tified in ref. [26J. Namely, (a) the superluminalilty of scalar excitations about some physically 
relevant solutions, and (b) the excessive softness of certain scattering amplitudes at low energies. 
These features are in fact related to each other, and although they do not necessarily signal any 
inconsistency of the low-energy effective theory — unless closed time-like curves may form — they 
do indicate that our model cannot arise as the low-energy description of a Lorentz-invariant mi- 
croscopic theory obeying the standard axioms of S-matrix theory, like e.g. any renormalizable 
relativistic QFT [26J. This is to be contrasted with what happens for GR. The Weinberg- Witten 
theorem makes the quest for a UV-completion of GR within local Lorentz-invariant QFTs particu- 
larly hard — perhaps hopeless [28] • Nevertheless, looking at GR as a perturbative Lorentz-invariant 
effective theory of gravitons living in Minkowski space, one finds that gravitational interactions 
cannot lead to superluminal propagation with respect to the underlying Minkowski light-cone as 
long as sources obey the null energy condition [291 [26] . Therefore, although probably not local, 
the UV completion of GR may well be Lorentz-invariant (like e.g. string theory). 

From a more theoretical perspective, we have constructed a class of models with peculiar 
derivative interactions. On the one hand, 7r's self-interactions involve many derivatives, more 
than one derivative per field. On the other hand, the particular combinations we isolated yield 
second-order equations of motion. At the classical level, this means that no new ghost-like degree 
of freedom shows up in the non-linear regime. This is a necessary condition for the quantum 
effective theory not to break down for sizable classical non-linearities. In fact with a consistent, 
if not generic choice of counter-terms, our effective theories can be extrapolated deep into the 
non-linear regime [20]. In a forthcoming publication, we show that these peculiarities allow our 
models to violate the null energy condition with no manifest pathologies in the low-energy theory. 
In particular, our models subtly evade the theorems of refs. [7J [32] . 
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Appendix 

A Existence and uniqueness of galilean- invariant terms at 
each order in tt 

As we discussed in section [T^t the generic structure of a galilean- invariant lagrangian term with n+ 
1 7r-factors is (d 2 TT) n ~ 1 diTdiT and the corresponding equation of motion has exactly two derivatives 
acting on every tt, (ddix) n = 0. However, because of the shift symmetry of the tt Lagrangian, 
the e.o.m. for tt can also be written as the divergence of the conserved current associated to 
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this symmetry, schematically: d fJ, (d 2n ~ 1 TT n ) fl = (dd7r) n = 0. This equivalence can be read as the 
statement that the combination (ddiT) n that appears in the e.o.m. derived from a galilean-invariant 
term is a total derivative. 

In the following we show that there exists one and only one total derivative (dd7c) n for every 
n; for this reason there can be at most one galilean-invariant term at each order and eventually 
we conclude that such a term exists for every n by constructing it explicitly, starting from the 
corresponding total derivative. 

First of all, notice that it's easy to write down a total derivative at the n-th order: 

. . . rT pK) . . . d, n d Un n = T^ -^d^n . . . d, n d Un 7i , (132) 

p 

here the sum runs over all the permutations p of the indices ^ and (— l) p is the parity of the 
permutation. We now prove that this is the only one for a given n. In what follows we use the 
same notation of section 11.41 II is the matrix of second derivatives of tt, = d^d u ir and the 
brackets [M] denote the trace of the martrix Ai, hence [IP] is the cyclic contraction of order n. 

Supppose that there are two total derivatives with the same n. Take the linear combination L n 
of the two such that the terms with the highest number m < n — 1 of Uhr are c^ l (Cl7r) m [II* 1 ] . . . [IP*] 
with i\ + . . . + ik = n — m and ij > 2 for j = 1, . . . , k: actually we can always cancel in C n at least 
the term (Chr)™. The multi-index % counts the different contractions that can appear in (ddTr) n ~ m . 

Now if we consider the combination L n as a Lagrangian density, by definition of total derivative 
the equation of motion for 7r that we get from C n must vanish. Let's see what are the consequences 
of this condition on the coefficients c l m . To find the e.o.m. we can integrate by parts twice and 
isolate a ir shifting the two derivatives acting on it to all the other II in turn, in particular when 
we choose a 7r in the j-th cyclic contraction and we let the two derivatives act on the 7r close to the 
one we are starting with in the same cycle, we end up with the structure [CM . . . IP 3 ~ 2 ]. Therefore 
among many terms, C n will give a contribution to the e.o.m. proportional to 

^(□TrniT' 1 ] . . . [ir*- 1 ] [an. . .n*- 2 ][n* +i ] . . . [u ik ]. (133) 

In order to have a vanishing equation of motion, this must be cancelled by a contribution coming 
from another term in C n . The only possible contraction that gives something proportional to 
(I133p in the equation of motion would be given by 

C n D (□7r) m+1 [n il ] . . . [IT^- 1 ] pT*^ _1 ] [IT*+ 1 ] . . . [IP fc ]. (134) 

In facy if we integrate by parts and we shift the □ of the (m+l)-th \3tt to one of the II in the j-th cy- 
cle we get exactly tt . . . f!133j) : this is the only other possibility to have in the e.o.m. a □ in the struc- 
ture [CM . . . IP J ' -2 ], maintaining at the same time (□7r) m and the cycles ii, . . . , ij-i, ij+i, ■ ■ ■ ,ik- 

However, terms with (□7r) m+1 like fjl34j) are absent from C n by construction, and as a conse- 
quence C n can be a total derivative only if c l m = 0. By induction on m we conclude that the whole 
C n must be zero: the two total derivatives we started with are the same a part from an irrelevant 
overall constant factor. In other words, the structure (ddn) n can be a total derivative only if it 
contains \3n to all powers from to n, and this can happen only once for a given n. 
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Since we proved that the only total derivative for every n is given by fll32j) . the last step of 
our argument is to show that this specific contraction can be always derived as the equation of 
motion of a galilean-invariant term. The galilean term can be constructed by replacing in turn 
each derivative with d^/iv: 

^galilean = T^-^d^d^ d^TT . . . d^TT + . . . + d^TT . . . d^_, d Vn _ x 7T 8, n 7T C^Tr) . 

(135) 

It's easy to check that the equation of motion we find starting from (11351) is precisely — 2n p^-^ 
d^d^iT . . . d^ n d Un 7T = 0; in fact terms with 3 or 4 derivatives acting on the same 7r must have 2 
derivatives with a v index like d Vi d Vj d^ k Ti or d Vi d v d ilh d il {K and because of the antisymmetry of 
rppxm—iinvn j n exchange of i/j <-> Vj due to the factor (— l) p , they are all zero. 

B Q <-> di 

We want to expand the equation of motion about the deSitter solution, 

n^-\H 2 + 7i , d^d v 7T -> -iffV + W** ■ ( 136 ) 

Each Galilean term (|3"9TT4"3"1) will then yield a linear combination of itself and the lower-order ones. 
The contributions to the tadpole E\ coming from all the terms will cancel — since we are assuming 
that the deSitter background is a solution to the e.o.m. We will therefore neglect the tadpole and 
concentrate on the other terms. Upon the replacement (11361) we get 



£ 2 -> £2 (137) 

£ 3 -> £ 3 -3H 2 £ 2 (138) 

£4 -> £4 - 3H 2 £ 3 + \H 4 £ 2 (139) 

£ 5 -> £ 5 - 2H 2 £ A + 3H 4 £ 3 - 3H 6 £ 2 (140) 



Therefore, comparing eq. (j4"5l) with eq. fj4"9~]) we have 
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The inverse relation is 
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(142) 



Notice that the two matrices that describe the mappings {c} — > {d} and {d} — > {c} are exactly 
the same, upon replacing if 2 with —if 2 . This had to be expected, for we can think of the trivial 
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configuration as an excitation of the deSitter one with 7T = + \H 2 x 2 . Then, like in ( 1T36D . to get 
the equation of motion for 7r fluctuations about the trivial configuration we have to perform the 
replacement 

7T -> +\H 2 X^ + 7T (143) 

everywhere in the n e.o.m. about the deSitter solution. Of course this way we will not recover the 
tadpole term, which we can however reconstruct by imposing the existence of a deSitter solution 
with given H 2 . 
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